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Abstract 

We address the question whether the cut-off dependence, which has to be introduced 
in order to properly define the Lippmann-Schwinger equation for the one pion exchange 
potential plus local (5-function) potentials, can be removed (up to inverse powers of 
it) by a suitable tuning of the various (bare) coupling constants. We prove that this 
is indeed so both for the spin singlet and for the spin triplet channels. However, the 
latter requires such a strong cut-off dependence of the coupling constant associated 
to the non-local term which breaks orbital angular momentum conservation, that the 
renormalized amplitude lacks from partial wave mixing. We argue that this is an 
indication that this term must be treated perturbatively 
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1. Introduction 

Since the original suggestion by Weinberg |j] that the nuclear forces could be understood 
within the framework of effective field theories (EFT) there has been an increasing interest 
in the subject (see for recent reviews). A key ingredient of the EFT formalism is that the 
cut-off dependence which is introduced in order to smooth out ultraviolet (UV) singularities 
can be absorbed by suitable counterterms, and hence any dependence on physical scales 
much higher than the ones of the problem at hand can be encoded in a few (unknown) 
constants. In order to achieve this in a systematic manner counting rules are also necessary. 

Weinberg's suggestion consisted of two steps. The first one was calculating the nucleon- 
nucleon (NN) potentials order by order in Chiral Perturbation Theory (%PT) from the 
Heavy Baryon Chiral Lagrangian (HB^L) ||. The second one introducing the potentials 
thus obtained in a Lippmann-Schwinger (LS) equation. There is no doubt that the first step 
can be carried out within an EFT framework: the renormalized NN potentials are known at 
leading, NL and NNL order |4|, [5], with even higher order corrections || and isospin breaking 
terms ]7| taken care of. The second step however is delicate. The potentials obtained in the 
first step are increasingly singular at short distances as we rise the order of xPT they are 
calculated. Hence the introduction of a regulator in the LS equation is compulsory. Since, 
even with the leading order (LO) potential, the LS equation can only be solved numerically, 
it is not clear that the scattering amplitude thus obtained is cut-off independent. This is so 
even for the successful fits ||, ^| to different partial amplitudes, where the cut-off is regarded 
as a variational parameter close to the last scale integrated out. We present here a proof 
that this cut-off can be removed from the LO (in the %PT counting) NN interaction if we 
tune properly the coupling constants of the potential. However, for this to be so we also have 
to tune the coupling constant of a non-local potential in the triplet channel. Even then, the 
only solution we find turns out to be physically unacceptable. Nevertheless, the insight on 
scaling so gained enables us to put forward a new proposal of counting rules where, coming 
back to standard procedures, divergences are fully absorbed by local counterterms. 

As EFTs have been mainly used in a perturbative framework, it is far from obvious 
how the two main features of those, namely renormalizability and counting rules, must be 
implemented in a non-perturbative one. Although in this work we shall primarily address 
the question of renormalizability, we would like to start by making a remark on counting 
rules, which emanates from previous experience on EFT in non-perturbative systems. It 
was pointed out in ref. || that calculating the potential in a non-relativistic system can 
be understood as the integration of certain degrees of freedom, which can be implemented 
as a matching calculation between two EFTs. In our case the higher energy EFT is the 
HB^L for the two nucleon sector, which is a local theory with pions and non-relativistic 



nucleons as explicit degrees of freedom. This EFT has an energy (E) cut-off (A#) such that 
E ~ m^ « A E « M ~ ill f n and a momentum (p) cut-off (A p ) such that p « A p « 
M ~ 4^/^ (m n and M stand for the pion and nucleon mass respectively, and 47r/ 7r for the 
scale of non-Goldstone boson QCD states). Its lagrangian can be organized according to the 
chiral counting since chiral symmetry (and its breaking) are explicit. The lower energy EFT 
has an energy cut-off (A' E ) such that E « h.' E « m n and a momentum cut-off such that 
p < m n << A^ << M . It consists of non-relativistic nucleons interacting through a (non- 
local) potential. The potential plays the role of a matching coefficient. As such, the potential 
encodes information on the higher energy EFT and it can be calculated independently on 
how the calculation of the lower energy EFT is organized, namely independently on what the 
counting in the low energy EFT is. Hence, on the one hand, the potential can be calculated 
order by order in %PT. On the other hand, chiral symmetry is not explicit anymore in the 
lower energy EFT (no pion fields exist) and consequently the chiral counting is not the 
natural way to organize the calculation anymore. 

An interesting example of a related situation is the pionium system (see |1 for a recent 



account), which has been studied using a series of EFTs [|10J- The higher energy EFT is 
the Chiral Lagrangian coupled to electromagnetism and the lower energy one a quantum 
mechanical Hamiltonian with the Coulomb potential and local interactions. The matching 
between the two EFTs can be carried out perturbatively in %PT and a, but the calculations in 
the lower EFT are carried out keeping the Coulomb potential non-perturbatively (otherwise 
no bound state exists) and, furthermore, one does not need to specify to which order of x?T 
the local potentials have been calculated. 

Following that spirit, the main question for the NN system is what should be treated 
as the LO potential in the low energy calculations. In the (higher energy) %PT counting 
the LO potential consists of the one pion exchange term (OPE) plus two local (^-function) 
terms. This assumes that the natural scale of the two local terms is of the order of the last 
scale integrated out (~ M). If the NN system was in a perturbative regime the scale of these 
two local terms would provide the scale of the scattering lengths. Since the experimental 
scattering lengths are much larger than the ones predicted in this way, we can foresee at least 
two possibilities. The first one is that an unsuspected behavior of QCD at energies ~ Aqcd 
produces unnaturally large values for the local potentials. Then one may consider these local 
terms as the (low energy) LO potential and treat the OPE (and higher orders) perturbatively 
nfl . This approach has been worked out at N 2 LO |I2] showing slow-convergence in the 1 5'o 
channel and no convergence at all in the 3 Si- 3 D\ channel. The second possibility is that the 
local terms do have natural sizes but the low-energy dynamics is responsible for the large 
scattering lengths. In this case there is no reason to treat the OPE perturbatively and a fully 
non-perturbative evaluation of the LS equation with LO potential (in the ;\TT counting) is 



required ||, ||. We shall stick to this second possibility for most of the paper, although 
eventually a third possibility, which is half way, will emerge as the most reasonable one (to 
us). 

Before going on, let us briefly discuss some previous work on the renormalization of the 
LS equation. The case of a local (and hence separable) potential, namely consisting of delta 
functions and its derivatives, has received plenty of attention |[|, [TJ], |I3| . This was expected 
to mimic the very low-energy (p << m^) behavior of NN scattering. The regularization 
of this pure local EFT was a matter of debate some time ago: a cut-off regularization 
showed a systematic order by order improvement in the phase shift fit whereas dimensional 
regularization (DR) with MS scheme was extremely sensible to the large scattering length 
and shallow (nearly)-bound state, which translated into a poor radius of convergence. The 
shortcomings of DR with MS were cured using the PDS scheme [nj (see also [[RJ). The final 
outcome appears to be equivalent to the well-known Effective Range Expansion [0. The 
next step in difficulty is renormalizing the LO potential in the 1 5'o channel, which contains a 



non-separable piece from the OPE. It was first carried out in [Q , and reproduced by several 
authors (see |19 , for a recent report). We shall re-obtain these results in section 3. Finally, 
as for renormalization in the 3 Si- 3 Di channel, the available literature is, on the contrary, 



somewhat scarce |19|, |2(J and the results are, to our understanding, not fully satisfactory 
(see sec. 6). 

The main difference of our approach with respect to the previous ones is that, in addition 
to the bare constants associated to local terms in the potential, we will also allow, but only in 
a initial stage, the bare constants of the non-local potentials to have non-trivial flows Q This 
is less restrictive than the standard assumption that only local terms should renormalize the 



LS equation [3^j for NN systems, which, in any case, is contained in it. This will permit us, 
not only to make meaningful comparisons with related work, but also to draw restrictions 
on the power counting. Having examined which conditions on the coupling constants are 
required in order to renormalize the LS equation, and the eventual consequences this has 
on our observables, we will be able to glean which terms of the potential can be included at 
LO, and which ones must be treated as perturbations, in the relevant case where only the 
coupling constants of local terms are allowed to flow and non-local potentials are fixed at 
the HB;\PT values. 

So, once the (low energy) LO potential has been identified, we suspect that, in order to be 
renormalizable, a higher order calculation should be organized as follows. The LS equation 



1 In fact, this turns out to be the usual approach in theoretical works on renormalization of singular poten- 
tials (see for instance |2l|]). Besides, there are known examples in a non-relativistic EFT of QCD (pNRQCD) 
where the renormalization of non-local potentials is required in order to absorb certain divergences |22[ , the 
most spectacular of which being the renormalization of the static potential [E3[ . 



must be solved and renormalized treating the LO potential, as well as its couplings, non- 
perturbatively, but the NLO potentials and higher perturbatively. 

Therefore, the first step in this program is to identify a (low energy) LO potential and 
to prove that it is renormalizable. We start by the naive choice, namely the LO potential in 
the xPT counting. For the spin singlet channel the LS equation is indeed renormalizable. 
However, for the spin triplet channels, we find that the LO potential in the x?T counting is 
only renormalizable if a certain coupling constant of a non-local potential has a non-trivial 
flow. (Or, in other words, if only the coupling constants of the local potentials are allowed 
to flow it is non-renormalizable.) Even in this case, the physical outcome is not satisfactory: 
the partial wave mixing is washed out of the renormalized amplitude. We conclude that the 
(low energy) LO potential is not the full LO potential in the x?T counting for the triplet 
channels. We identify a (low energy) LO potential, which is renormalizable, and prove that, 
if we treat the difference as a perturbation, it is also renormalizable at first order. 

We distribute the paper as follows. In section 2 we introduce a convenient basis for the 
NN wave functions and our notations. A brief note at the end of this section serves to close all 
what refers to the isosinglet-singlet channel. In section 3 we prove that the isovector-singlet 
channel is renormalizable and provide explicit expressions for the cut-off dependence of the 
bare parameters both for a hard cut-off and for dimensional regularization. In section 4 we 
prove that the isosinglet-triplet channel is also renormalizable, but requires a strong cut-off 
dependence of the coupling constant of the (non-local) term in the potential which, in turn, 
prevents the renormalized amplitude from partial wave mixing. We interprete this result as 
an indication that this term must be treated perturbatively and prove that, if so, the first 
order in perturbation theory is finite. After briefly discussing in section 5 the isovector-spin 
vector channel, section 6 is devoted to a discussion. 

2. A convenient decomposition 

We start from the LO NN potential given for instance in ref. ||: 

T//1, l^ f 9A \\ ^i ■ (k - k') a 2 • (k - k') 

V(k,k) = -[jj-J n-r 2 {k _ k , )2 + ml + C s + C T a 1 .a 2 , (2.1) 

where f n is the pion decay constant (~ 93 MeV). 

This potential acts on a wave function \l/^(k, k'), where a, b and a, (3 are nucleon isospin 
and spin indices respectively. This wave function can be decomposed into irreducible repre- 
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sentations of spin and isospin as follows: 

\r 2 ) ab (a 2 ) a ^ ss (k) + (T 2 ) ab (a k ,a 2 ) a ^ k s V (k) + 
+(r k T 2 ) ab (<T 2 U^ s {k) + (T k T 2 y b (a k ,a 2 ) af3 ^ k v k ;(k) 

The potential Q2.1|) reduces for each isospin-spin channel to: 
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We still have to implement Fermi symmetry. This implies that the irreducible wave 
functions ( |2.2|) must fulfill (isospin and spin indices will be omitted for the rest of this 
section): 

ipss(k) = -ipss(-k) , 

^sv(k) = ^(-k) , 

Vvs(k) = ip vs (-k) , 

Vvv(k) = -^yy(-k), 

which is implemented in the LS equation if we choose: 

T ss (k, k'; E)= l - (V ss (k, k>) - V ss (-k, k')) + 
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T sv (k, k'; E) = - (V sv (k, k') + V^(-k, k')) + 



1 /" A rf 3 fc" 
+ 2~J (2^ 
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It is the advantage of the above decomposition that we will not need to specify which 
(coupled) partial waves we are analyzing. 



If the LS equation for the potentials fl2.3|) was well defined, using ( |2.5| ) would be equivalent 
to solving the LS equation: 



T xy (k, k'; £) = V xy (k, k') + / -—3 \4,(k, k") — -^ T xy (k", k'; £) , (2.6) 

(x,y=S, V), namely ignoring the statistics and then using the standard formulas: 

T ss (k,k';E) = ^(f ss (k,k';E)-f ss (-k,k';E^ (SS^VV), 

T sv (k,k';E) = ^{f sv (k,k';E) + f sv (-k,k';E^ (SV^VS). (2.7) 

However, the LS equation for T xy is not well defined in any channel and hence using (|2.5| ) 
or fl2.6p-(P77p may not be totally equivalent. In particular, for the SS and VV channels, the 
UV divergences one finds using ( |2.5| ) are softer than those from fl2.6|) -( p^TD , so we shall work 
with (|2.5|) . For the SV and VS channels, however, the UV divergences found using (|2.5|) are 
as strong as the ones that stem from ( [2.6| )-( f2~7] ). For convenience, we have chosen to work 
with the latter for these channels. 

The LS equation in the isoscalar-scalar channel in (|2.5| ) is already well defined, as it is 
apparent from the anti-symmetrization of the corresponding potential Q2.3p . On the contrary, 
the other three channels require regularization. Searching for the systematics to tackle them 
will be the aim of the next three sections. For notation simplicity, the energy dependence 
of the T-matrices as well as of other auxiliary functions will not be displayed explicitely for 
the rest of the paper. 

3. The isovector-singlet channel 

The LS equation for this channel reads: 

6 

(2*) 
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where in the last lines we remember the values those constants would take if the potential 
had been calculated at LO in %PT. The hat and the VS subscript will be dropped in the 
following. 

Let us define: 



A(k')8 ij := I 
Then Q3.1|) reads: 
T ij (k, k') = co(l+A(k'))8 ij + 
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(3.4) 



T(k,k') = c 2 T 2 (k,k') + c (l + .4(k')) 
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where we have dropped the b %i structure. If A(k') was a fixed function, the equation above 
would be well defined and could already be solved with no need to regularize it. However 
•A(k') is a functional of T and a second equation which relates them must be introduced. 
This is achieved by multiplying eq.(|3.5|) by 1/(E — jt + if]) and integrating over k. We 
obtain: 
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Substituting iteratively T 2 in ( |3.4j ) in the rhs of ( |3.6| ) we see that only the first iteration 
produces further divergent expressions when A —>■ oo. We can then write ( |3.6| ) as: 
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where Xq and C contain linearly and logarithmically divergent terms respectively, whereas 
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It is clear that the expression (|3.7p can be renormalized by a redefinition of cq. In 
dimensional regularization, (D=3+2e), we obtain: 

I M2 °2 fl i ^ 1 

Xms = 7£-Log(4vr), (3.9) 



which is in agreement with [18], and for a hard cut-off: 



1 MA M 2 c 2 /A 2 \ 

2 Log - +-— r. (3.10) 



c 2tt 2 32vr 2 \fJ? J c r (fi) ' 

If we now wish to solve numerically the LS equation, we should proceed as usual and 
introduce a hard cut-off. However Co is not to be fitted to the experimental data but substi- 
tuted by ( |3.10| ) and the cut-off made as large as possible (in practice it should be enough if 
VEM/A is of the order of neglected subleading contributions from the NLO potential, (see 
for a more technical discussion). What we have just proved is that the result will be 



cut-off independent up to corrections yEM/K.. fi must be fixed at the relevant momentum 
scale fi ~ (vEM,m n ) and Cq(/x) tuned to fit the experimental data. 

Although we have no prediction for c^lfi) we can try to understand from ( |3.10p how 
large scattering lengths may arise. Since Cq(/x) evolves according to a non-perturbative 
renormalization group (RG) equation it might take very different values depending on the 
scale it is evaluated at. After solving it: 

d&p) = -^M . (3.11) 

i + M y Lo gA 

1671-^ ° /jo 

if we input the value of ref. |L8| c^m^) = — ( 79 ^ eV ) 2 , we obtain Cq(M) = — ( 125 m gV ) 2 ? which 
is not quite at the natural scale (~ M). Hence, the non-perturbative low energy dynamics 
does not seem to be enough to fill the gap between the natural scales and the large scattering 
lengths. In spite of that, the variation of Cg(/x) from m w to M is large enough as to justify a 
non-perturbative treatment of the OPE in this channel. 



4. The isosinglet-vector channel 

The LS equation for this channel reads: 
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where we show also the LO values of the coupling constants. We shall drop the subscript 
SV and the hat in the following. We call the term proportional to c\ above spin symmetry 
breaking (SSB) term. This term breaks orbital angular momentum conservation and makes 
the analysis of this channel qualitatively different from the previous one. In order to illustrate 
it, let us take k' = for simplicity. As we regulate (|4.1|), the possible divergences arising when 
the regulator is removed depend on the high momentum behavior of T y '(k). If T iJ (k) ~ |k| a , 
the usual power counting arguments imply that, due to the SSB term, the integral on the rhs 
will rise this power by one. Hence, the high momentum behavior of the lhs of the equation 
will not match the one of its rhs unless: (i) a = — 1 and the high momentum contribution of 
the potential cancels out the one arising from the integral or (ii) a = and the bare coupling 
constant c\ goes to zero as the cut-off goes to infinity, which removes the |k| Q+1 term on the 
rhs. We prove in the Appendix A that the case (i) in fact reduces to (ii). 

The preceding discussion provides a rather intuitive introduction to what, in the course 
of section 4.1, we will demonstrate in full detail. That is, all those rising divergences caused 
by the SSB term can only be renormalized by a, so far undetermined, flowing of their 
accompanying coupling constant, c\. Next we will fix this cut-off dependence and, having 
explored the consequences such a behavior has on the amplitude, will come back in section 
4.2 to standard procedures. There it is shown that the alternative of treating SSB as a 
perturbation solves the problem, as all divergences get renormalized by local counterterms 
and no C\ flowing is longer required. 



4.1 Non-perturbative treatment of the SSB term 

Let us then return to equation ( |4.1|) . It has the following structure: 

(k-k'Wk-k'V 



r j (k,k') = co(^ + ^(k')) + Cl 

5 ij 



(k - k') 2 + m 2 



+ & 3 {k,k') 



+ 



+C 2 



+ C 2 



A ij (k') 
i^(k,k') 

Let us define: 



(k-k') 2 + m2 

A d 3 k T^'(k,k') 
^YE-^ + i-q' 

A J3iU// 



d 3 k'' 



T l Hk\k' 



(2tt)3 (k - k") 2 + mlE-^ + 



M 



17] 



d 3 k" (k-k") i (k-k") fc T^(k",k') 
(2vr)3 (k-k") 2 + m2 E -^ + 



M 



IT] 
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which allows us to isolate in T-^(k, k') and Co (5 l: > + A 1 -* (k!)) all sources of divergent behavior, 
since T (k) and T 2 (k, k') are perfectly well defined. 

Using the expressions of B ij (k,k') in (|4Xl| ) and T ij '(k,k') in (gXgD , Tf(k,k') can be 
re-casted in the form: 
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This decomposition enables us to compute cq(5^ +*4 y (k')), and hence the full amplitude 
T lJ (k, k'), in terms of T (k), T^(k, k') (n — 0, 1, 2) and T 2 (k, k') through the equation 
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As we have already mentioned To(k) and T 2 (k, k') are finite when the cut-off is removed. 
If we solve T^(k, k'), n — 0, 1,2 iteratively, the most divergent pieces in the n-th iteration 
are Tio ~ (ciA) n A, Tn ~ (ciA) n and T i2 ~ (ciA) n_1 ci. These series are expected to have a 
finite radius of convergence. The radius of convergence is in any case non-zero because they 
are bounded by geometric series (or derivatives of them). If c\ does not go to zero as 1/A 
or stronger (in particular, if C\ is not allowed to flow), each series will separately diverge. 
In that case, a finite result can only be obtained if non-trivial cancellations occur for all n, 
which we do not see how they could actually happen. If, on the contrary, 

Cl (A) = °j- + .... (4.1.5) 

and C\ is small enough, the series will converge. For the T-matrix, such a strong cut-off 
dependence implies that the terms: 

Cl T^(k,k') — *£2^+..., 

Cl rr 2 (k,k') — f M^l+..., (4.1.6) 

where t 10 is simply a finite constant and, as we see, all k, k'-encoded information will be 
washed out from the amplitude. 
That is to say: 

T\k, k') = Urn c (S kj + A kj (k>)) (T (k)cf fc + Cl Ti k (k)) + Cl Tl{(k, k') + c 2 (T 2 (k, k') 8 ij + 
+ Cl TH(k, k')) = c (S ij + A l >{k!)) (r (k) + Cl tg) + c 2 T 2 (k, k')^' , (4.1.7) 

which is finite provided Co (5^ + A^ik!)) is finite. In order to prove the latter we borrow 
from section 3 the following results: 

A r„(k) _ MA + A^ Log ^ +0(l) _ 



(2i) 3 E-f + ii) 2tt 2 32tt 2 b Va* 

A d 3 k T 2 (k,k') 



(2tt)3 £ _ g + 



a; 



277 



0(1), (4.1.8) 
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and find in Appendix B: 



Cl / (t t F T i (k) ■ = aoA + ifooV^MH-doLogfAUo^ 

y (2tt) 3 £ - f- + i V Xm^J \A 



(2tt) 3 


£- 


17 + ^ 


A d 3 k 


rpil 

1 n 


(k, k') 


(2tt) 3 


E- 


k 2 . ■ 


A rf 3 A; 


rpii 
1 12 


(k, k') 



ci ,„ ' _ "v = o(i: 



where a , &cb ^o are cut-off independent constants related to c.\. Then the flow: 

1 MA,a„A MV M 2 \ rfo T M 2 \ 1 r/l 1in , 

+ ^zrLog -r + -Log — + — -r , (4.1.10) 



c 2tt 2 3 32tt 2 \fi 2 J 6 V/x 2 / c^(/x) 



makes Cq{5 1 ^ + ^.^(k')) finite and hence does ( 4.1.7| ). We have then proved that the flows 



( [4.1.5 ) and (|4.1.10| ) renormalize the triplet channel. 



It is not difficult to see that the various series above involving divergent terms are bounded 
by geometric series or derivatives of them. This ensures that our flows provide actually finite 
expressions for the amplitude if C\ is small enough. However, this amplitude appears to be 
diagonal in spin space and hence orbital angular momentum is conserved. Although, the 
observed 3 Si- 3 Di mixing, which is small, might be attributed to a higher order effect, it is 



clear from ref. [25| that it is due to the OPE to a large extend. In order to preclude the 
conservation of orbital angular momentum, we can foresee two ways out: (i) a SSB term 
may survive in the renormalized amplitude if c.\ is tuned infinitely close to the radius of 
convergence of the series, so that our bounds do not hold anymore, and (ii) the SSB term 
from OPE must be treated as a perturbation and renormalized as such. The possibility 
(i) is examined in Appendix C, where we show it unlikely to be realized. In the following 
subsection we explore (ii) and prove that if a suitable SSB term is treated as a perturbation, 
the amplitude is renormalizable at first order and the mixing survives. 
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4.2 Treating the SSB term perturbatively 

Let us split the potential as: 

V ll (k,k') = V {0)ij (k,k') + V {1)ij (k,k' 
V^ )l Hk,k') ic + ^ r4 At ij 



(k-k') 2 + m2 

k-k'Y(k-k'y -^^5^ 



2 

A2 



with LO values for the coupling constants indicated. In the following we drop the SV-channel 
sub- indexes. 

The amplitude will be written as: 

T ij (k,k') = T {0)lj (k,k') + T (1)ij {k,k'), (4.2.2) 

where T(°)^'(k,k') fulfills: 

/A ,73 L// i 
i^_0o)ifc( k;k ") i T^ kj (k",k') . (4.2.3) 
( 27 E-\r+iri 



M 



The renormalized solution to this equation is given by T^ y (k, k') = T(k, k')5 lJ ' in section 
3. At first order in perturbation theory r^ tJ '(k, k') verifies: 

T^'(k,k') = V^ ij (k,k') + ^y« ifc (k,k") -4 T(°)^(k",k') + 

J ( 27r ) E-^ + ir] 

^^k") - JL . TW(k»,k>) . (4.2.4) 

( 27r ) E-^ + ir] 



Using (|3.4j ) and (|3.5|) we can see that the second term above is finite. We can then gather 



the first and second terms into a new, energy dependent, potential defined as: 

~ r A rl 3 k" 1 

V^ ij (k,k") :=^ (1)ij (k,k')+ / j^V {1)ik (k,k") ~ T^ kj (k",k'). (4.2.5) 

J ( 27r ) E-^ + ir] 
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Therefore, the integral equation reduces to: 

d 3 k" Zo twhwm 



t^M) = ^(k,k'0 + ^(k + / {2nr(k _ kr + mlE _^ + 



TZ ij (k') := J 



M ' V 

A 



d 3 k » T (i)y( k /' )k /) 



(4.2.6) 



In order to prove it finite we decompose: 



T (1)iJ (k,k') = coTZ kj (k')T* k (k)+Ti J (k,k'), (4.2.7) 

with T y (k) defined in ( 4.1.2Q and T[ 3 (k, k') given by: 



~ ~ f A rl 3 k" r„ T ij (]<:" y\ 

T^X> := ^( k , k ») + / ^-^JL|^ . (4 ,,) 



Both T y (k) and T^ (k, k') are well defined (the tensor structure is crucial for the latter 
to be so). Divergences can only arise in coT^^k'), which reads: 






J (2tt) 3 p k- 

The numerator is well defined (for that the tensor structure is again crucial) and the 
divergences in the denominator have exactly the same structure as in the denominator of 
( p.7|) . Hence they are renormalized by the same cq flows. We have then proved that if 
we treat the SSB term as a perturbation, the amplitude is renormalizable at first order in 
perturbation theory and no extra counterterm needs to be introduced. 

5. Isovect or- vector channel 

If we use (|2.6|)-( |2~7D in order to obtain Tyv(k, k'), the calculation of Tyy(k, k') would 
reduce to that of the previous section. However, as mentioned in section 2, the UV behavior 
is smoother in terms of (|2.5| ), as it happens in the SS channel, although here we still need 
to introduce a regularization. The LS equation, dropping de isospin delta, reads: 

"A rl3t,ll 1 

M+ ir ) 

where: 

v A '-(kk>) - 1 (v- (k k') v- ( k id)) - Cl ( {k - k ' nk - k ' y (k + k ' Y{k + k ' y 

Vyy (k,k) - 2 \V V y{k^) V vv { k,k)j- 2 ^ (k _ k/)2 + m 2 (k + k0 2 + m 2 

I c > ( 5 1 5 1 \ (5 2) 

2 V(k-k') 2 + ^' {k + k'f + ml) ' K ' ' 
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/A j3 U» I 

-— ¥ l#*(k, k") k „ 2 ■ T&(k", k') , (5.1) 

( 27r ) E - \r + in 



where those constants calculated at first order in x?T take the values: 

We have not analyzed the possible existence of non-trivial flows which may renormalize 
the above equation. The fact that the SSB term must be treated perturbatively in the SV 
channel, indicates that also here we should proceed according to the same philosophy. The 
potential ( |5.2| ) in the zeroth order approximation reads: 

which leads to a well defined LS equation. At first order in perturbation theory we will have: 

T^(k,k') = T^'(k,k')+T^(k,k'), 

/A j3i,h rp(0)kj /,„ , ,x 

+r^^"(k,k") ^if' k,) , (5.5) 



which is also well defined. We expect the divergences arising at higher orders to be absorbed 
by local counterterms. 

6. Discussion 

We have addressed the renormalization of the LS equation for the LO potentials (in the 
xPT counting) of the NN system in all channels. In addition, for each channel we have been 
able to carry out our analysis for all partial waves (including partial wave mixing) at once. 
The isoscalar scalar channel does not require regularization. For the isovector scalar channel 



we recover the flows of ref. [IS|. The remaining two channels have deserved a more detailed 
study. 

The first non-trivial result is that the renormalization of the isoscalar vector channel 
requires a strong flow of the coupling constant of a non-local potential, the SSB one, or, in 
other words, if only the coupling constants of the local potentials are allowed to flow, the 
isoscalar vector channel is not renormalizable. Several comments are in order. 



First of all, the flow ((4.1.5[) of the coupling constant of the SSB term is not such a big 



surprise. Notice that at high momentum this term tends to a (direction dependent) constant, 
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which is the same behavior (except for the direction dependence) as the 5-function term both 
in the singlet and the triplet channel, the coupling constants of which also show similar flows. 
The main difference is that the leading behavior for cq is fixed and the subleading contains 
the free parameter (eg (//)). For c\ instead, the leading behavior contains the free parameter 
(ci) and the subleading behavior is not observable. 

What is worse, the flow ( |4.1.5| ) has undesirable consequences: the renormalized T-matrix 



conserves orbital angular momentum, even if the bare interaction does not (see Appendix C) 
0. Since it is precisely the OPE the main responsible for mixing (also of higher partial waves 
p5|), we would like it to keep doing this job for us. Therefore, in order for c\ not to flow, but 
to be fixed at the HB%L values and produce partial wave mixing, we are forced to exclude 
the SSB term from the (low energy) LO potential, and to treat it as a perturbation. This 
also appears to be reasonable from the phenomenological point of view since the observed 
mixings are small [p5fl . 

We have developed this line in sections 4.2 and 5. We have proved that at first order the 
vector channels remain renormalizable (at zeroth order the problem reduces to the one in the 
singlet channels, which are renormalizable). The picture which emerges is half way between 
| ll| , where the pions are treated perturbatively, and ||, [17| where the whole potential is 



treated non-perturbatively. The (low energy) LO potential is the part of the LO potential in 
the ;\TT counting which conserves orbital angular momentum. We are tempted to propose 
the following counting. The 0(Q n ) (n = 0, 1, ...) contribution to the NN potential in the x?T 
counting must be divided into two pieces: the one which conserves orbital angular momentum 
(SS) and the one which does not (SSB). The SSB terms keep their xPT counting but the SS 
ones are enhanced and must be counted as 0{Q n ^ 1 ). Only the LO potential OiQ' 1 ) must 
be treated (and renormalized) non-perturbatively. We have seen here that this proposal is 
theoretically consistent at next to leading order, and, in addition, it does not require any 
coupling constant of a non-local potential to flow anymore. It remains to be seen if it is still 
so beyond that order and, of course, whether it is phenomenologically successful. 

Let us finally comment on recent work on the subject [|1^, p0| . The authors in both 
references try to renormalize the triplet channel by adjusting the coupling constant of the 
5-potential only. Hence, according to our results both works should show a remnant cut-off 
dependence when the cut-off is large enough. Note also that it is only in the large cut-off limit 
when a meaningful comparison is possible, since the regularizations used in the three works 
are different. The authors of ref. PC| , who use a subtracted (/i-dependent) LS equation, 

2 We have also checked perturbatively in c\ and ci up to order c\C2 that the effective range depends on 
c.\ only through the scattering length. Since the latter can be adjusted by tuning Cq(/i), up to this order 
both the scattering length and the effective range are blind to c\. We have not looked at what happens to 
the rest of the amplitude or to higher orders but we suspect that they are also insensitive to c\ . 
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argue that a reasonable boundary condition is that for large \x the T-matrix coincides with the 
potential, and check numerically whether, once the scattering lengths are fixed, the remaining 
observables are independent of \x for large ji. They find that for laboratory energies up to 
100 MeV. the 3 S\ and 3 D\ phase shifts are remarkably independent of /a for jjl > 0.8 GeV. 
but the mixing angle shows a strong /^-dependence for 6 GeV. > /a > 0.8 GeV. and only 
for fi > 6 GeV. the /i-dependence smooths and the results may appear to converge. We 
interpret this stronger /^-dependence of the mixing angle as an indication of the remnant 



cut-off dependence mentioned above. The authors of ref. [[19] obtain the flows by analyzing 
the short distance behavior of the Schrodinger equation (see also f26|| ). For the 1 S'o they are in 
qualitative agreement with ours. For the triplet channel they present analytic flow equations 
which are argued to coincide with those of the chiral limit. The flow of the S-i unction term 
is given implicitly by their equation (18). They assume that their q t , which is proportional 
to our Ci, does not flow^ and find a multi-branch structure for the flow of their Vo-R 3 , which 
is proportional to our cq (R — > 0, R playing the role of an inverse cut-off). It is interesting to 
note that if they allowed a n flow like our c\ in section 4.1, namely a n ~ R, and Vo-R 3 like our 
Co, namely Vo ~ 1/-R 2 , their eq. (18) becomes cut-off independent. Hence our flow (4.1.5) is 
a solution in the R — > limit to the flow equation (18) of ref. [[HJ. Recall, however, that, if 
oe n is not allowed to flow, the strict limit R — > cannot be taken. This is proved in Appendix 
D. Hence eq. (18) of Jl9[ does not produce an acceptable flow for V and, therefore, it cannot 
be used to properly renormalize the triplet channel. Furthermore, we would like to remark 
that leaving the cut-off finite and checking that the finite cut-off effects are higher order in 
the EFT expansion is a procedure which unavoidably will lead to problems in this case. If 
we wish to improve the accuracy of our EFT calculation, we will have to calculate at higher 
orders. Even if we insist in keeping R finite, we will have to choose it smaller and smaller for 
the LO terms not to jeopardize the accuracy of the higher order calculation. Then at some 
point R will hit the region were no continuous solution exists and we will loose all predictive 
power (if we give up continuity, an infinite number of inequivalent solutions exists). Note 
that the fact that finite cut-off effects can be compensated by higher dimensional operators 



27fl , which holds in perturbatively renormalizable (and asymptotically free) theories, needs 



not hold here. 
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A. The case a = — 1 

The more general decomposition of the high momentum behavior of T* J (k) for a = — 1 
reads: 



T ij {k) = B- 1 — + B_ 1 -— + V ij (k) 

k k r 



(a.i; 



where limk^ 00 'P u (k) ~ J%. Notice then that the integral in fl4.1|) at most diverges logarith- 
mically and, furthermore, the divergent term must be proportional to the <5 iJ tensor. By 
calculating the high energy behavior of the integral in the rhs of ( |4.1|) we obtain: 



T ij (k) 



c S ij + d 



W Mc 



o 



k 2 



Log(^)+A 



2tt 2 



B-! 



B-! 



Log 



A s 



-EM 



8 ij + 



Mci B-t + B-t 

4^" 3 



Mci 



£-i + 



B_i 



Log 



-EM 



+ / 2 






(A.2) 



with /i and ji two finite, constant terms. Observe that, although the cut-off dependence 
can be removed by a suitable redefinition of Cq, the non-analytic terms ~ Log|k| cannot be 
compensated by the potential. Self-consistency of (|A. 1|) and ( |A.2|) force c\ — > again. 

B. Proof of ( grt^l ) 



Let us define: 



H a (EM) := Cl /" 



ci 3 fc T&(k) 
(2tt)3 E _ g + irj 



a = 0,1, 2. 



and concentrate on Tio(EM) (the analysis for Tli(EM) is identical). We have: 



Ho(EM) 



d 3 k 



A A'iUH 



d 3 k" 



("1 



(k - k")*(k - k") j T (k") 



(2tt) 3 J (27r)^-g + ^ (k-k") 2 + m 2 E-^ + ir, 
x ,f'k f k d 3 k" Cl ci(k-k ff )*(k-k") fc + c 2 (J i * T^(k") 



(27T) 



(2vr) 3 £-£ 



M ' "7 



(k - k") 2 + ml 



E-^ + i V 



(b.i; 



(B.2) 



If we solve the equation above iteratively using ( |4.1.2j ) for T (k) and ( |4.1.3[ ) for T^(k), 
the most divergent term in the n-th. iteration is (super-indexes j and k are contracted with 



unwritten momenta): 

'n-i 

n 



r n+l 



H+2 M d 3 h 



(kx - k 2 ) i (k 1 - k 2 y 



1 



(2ir)* } E _ £ + ir} ( kl -k 2 ) 2 + m 2 

(k„+i — k n+2 ) (k n+1 — k n+2 ) 4 



E-§+i V 



1 



E-^ + irj 



»-n+l 



Hi+2 



,2 + m 2 



k 2 



(B-3) 



M ' "'I " "^ A/ + Z7 7 

Taking into account that the limits i? — > and m£ — > exist and the flow ( |4.1.5| ), the 
leading behavior in A reads: 

" +2 rA <i 3 A;/ 1 ' n ' - !- -^i' 1 ' - '' -Vj ^ - '' -^ 

/ ' 

Z=l 



M )n+2 c n+l J J-j 



(ki - k 2 )*(k 1 - k 2 )' 



(2vr)3k? 



(k 



k,) 2 



(k n +i — k n+2 ) (k n+1 — k n 
(k n+ i - k n+2 ) 2 



+2) 



~ Cl 



- X A 



which proves that ao is a (ci-dependent) constant. Notice also that the integral in 
is bound by (J <i 3 k/k 2 ) n+2 . Let us next identify the subleading behavior. Consider first 
E = 0. The derivative of ( |B.3| ) with respect to m 2 at m 2 = has at most a logarithmic 
singularity which means that the next to leading behavior in A is ~ c™ +1 A n_1 ?Ti 2 LogA, which 
gives rise to O (j^j contributions in (|4.1.9 ). Terms contributing to do in the nth. iteration 
appear when: (i) the c 2 -proportional term of T (k) is iterated through only C\ potential 
insertions coming from the second line in ( |B.2|) ; (ii) the equal to 1 term of To(k) is iterated 
in such a way that a c 2 potential from the last piece appears only once in the iteration. The 
relevant integral is obtained by substituting: 



(K 



kp+i 



(K 



V1) 3 



c 2 8^ 



(k p - k p+1 ) 2 + m 2 (k p - k p+1 ) 2 + m 2 

in (|B.3|) . In order to get the leading behavior in A of this integral we can set m 2 = in all 
but the substituted term above. We have (super-indexes j, I, q and k are contracted with 
unwritten momenta): 



A 



-my+ 2 c 2 c\{ n 

l=l\{p,p+l} 
A d 3 k p+1 (k p _i - k p )' 1 (kp_i - k 



d% 1 (k x - ka^kx - k 2 ) 



(2vr)3 k 2 



(ki 



k,) 2 






( kp+i — k p+2 ) 



1 (kp + i - kp 



(2vr)3 (kp_! 



kp) k p 



VH; 



+ m 2 



k P +i ( k P+i 



V2) 2 



c 2 a" LogA 



(B.6) 



which proves (with the flow (|4.1.1|) ) that d is a constant. 

Let us next address the energy dependent contribution to ( |B.1| ) . Notice that any analytic 
contribution in EM would show up at 0(1/ A). Hence only non-analytic contributions (like 



;b.4) 
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the one in ( P-6Q ) are relevant to us. Let us then look for non- analytic contributions in EM 
in the most divergent diagram in the nth. iteration (|B.3| ). Since the m% — > limit exists we 
can take it and have: 



n+1 f A dh kg f A dk n+2 kj +2 

. /\ n (ki - k 2 )'(k! - k 2 )* (k n+1 -k w+2 ) fc (k n+1 -k n+2 )* 
J (ki - k 2 ) 2 (k n+ i - k n+2 ) 2 

where dVti , i = 1, ..., n + 2 stand for angular integrals. Since the most singular contribution 
comes from the region |k;| ~ A VZ, the angular integral will give rise to a constant (which, 
furthermore, is bound by (47r) n+2 ), and the integrals over |k;| decouple. Hence the leading 
behavior for small E turns out to be the non-analytic contribution we are looking for («o, 
Po, 5o and (3q are constants): 

n+2 



„n+l 



A 

dk 



E-^+irj 



c™ +1 (a K + ifoVEM + O Q J j ~ 
^c^ +1 fa o A + o VE^ + o(j\\ , (B.8) 



which proves, in addition, that bo is a constant. Notice that a LogA dependence in this term 
would have been fatal for renormalization. 

We have then proved the first formula in Q4.1.9 ). The proof of the second formula is 



identical. The third formula is proved by simply noticing that all integrals involved are at 
most logarithmically divergent and, those which actually are, go multiplied by c\ ~ -?. 

C. On c\ tuning 

In section 4, when we focused on proving that a certain behavior of the bare constants of 
the potential as functions of the cut-off (namely, cq, c\ ~ A -1 ) would render a finite T-matrix, 
only c was conveniently fine-tuned. As a result the so-computed scattering amplitude 
lacked from partial wave mixing, which is expected due to the second rang tensorial term 
in the (bare) Hamiltonian. In order to obtain partial wave mixing two possibilities must be 
regarded. On the one hand, it could well happen that, indeed, mixing should not have been 
considered as LO, but as a NLO term to be treated perturbatively, the divergences it may 
cause being absorbed in the usual way by higher order local counterterms. This appears 
to be consistent with the fact that partial wave mixing in this channel amounts only to a 
few degrees. This treatment resums the (^-proportional part of OPE. Its SSB term, now 
eliminated by the strong suppression of Ci, is then recovered in a NLO analysis. We have 
shown how this works in the section 4.2 . 
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Nevertheless, another possibility remains unexamined. A proper tuning of c\ to a, let's 
say, non-trivial RG fixed point, could very well recover mixing at the leading order. So far, 
the existence of such a fixed point is anything but evident. Uncovering it or ruling it out 
requires detailed numerical work which is beyond the scope of this paper. However, in order 
to illustrate our point let us provide two approximations that exemplify how this tuning 
would emerge, how it would affect previous results and to which extent to achieve this goal 
we depend on the exact resolution of our actual system of integral equations. 

Let's take in the following k' = for simplicity. We will also apply the chiral limit (m^, 
C2 - > 0) and work with cq and c\ defined in section 4.2 . After decomposing the T-matrix in: 



T ij (k) = T^k)^ 



k l kP - ¥-6* 

k~2 



T 2 (k), 



(CI) 



the following two angular integrals arise in the resolution of its LS equation: 
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nk^nj _ h^Lfikj 
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'■^i F J k ^ L + Cl " 3 F 15 



(C.2) 



with uji (p-), i = 1, 2, 3, as known functions (k = |k|, k" = |k"|). 

At this point we wish to introduce some reasonable approximation that allows us to 
transform the non-separable in k and k" functions uii (p-) into separable ones. Once this 
is achieved we only need to solve a conventional system of equations and check whether, at 
least within this approximation, a non-trivial fixed point exists. Obviously our approximation 
should be as compatible as possible with what we know about the behavior of the full d 3 k- 
integrals. For instance: 



dk" k 



hi 



LU, 
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WE-fg + in 



k 



1,2 



(C.3) 



that is, both are finite integrals proportional to k in the limit A — ► oo. Unfortunately, no 
separable u>i achieves this. We shall content ourselves with a simple but still reasonable 
starting point that enforces separability. Then, let us take cj 3 (p-) as a constant (:= a^) 
and substitute uip (p-) by := a lj2 A (0:1,2 also being constants). Although the latter intro- 
duces logarithmic divergences which do not exist in the actual function, it keeps the correct 
behavior in k shown in (C3). 
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The LS equation takes the form: 

T*>'(k) = 00(1 + ^X0)^ + 0! 

k 2 



+cia 2 



k 4 k^ - fSv 



k 2 " 



+ Ciq;i 



k 2 



fcT 1 Z_ 1 + 



(C.4) 



where we have already used that Ti(A;) becomes momentum independent, as it is easily 
verified through ( |C.4| ): 



Ti = c (l + TiJ ) +Cia 3 i3, 
T 2 (/c) = ci(l + «i kTil-i + a 2 A;C) . 

The following functions have been introduced: 

; A dA; k 2 

X 

X-i 
C 



(C.5) 
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dk 




fc 2 T 2 (fc) 




2tt 2 
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c/A' 




k 




2tt 2 
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k 2 , ■ ' 


A 


(i/c 




kT 2 (k) 



A few manipulations allow us to solve for T\ and the combination ot\ T\X_\ + a 2 C: 

c + c^aaJo + A a i a z i-£ali 



(C.6) 



«i T{I_i + a 2 C 
where a quadratic divergence: 



iii-i 



1 - c X - cfaia 3 Tz^ a2lo 
ci« 2 X_i \ (c + cfoi 3 lo) ^ + ci(l - c X ) 



1 - Cl a 2 J ; 1 - Co J - 5 2 aia3 _Ai-_ 



(C.7) 



J i: = 



c/A' 



A- 8 



27r^-^ + ^ 



(C.S 



enters. 



It is not difficult to realize that little has been gained: the only way to get ( |C.7|) finite 
is by an untuned c 1; (1 — Ci« 2 X 7^ 0), and a tuned c , which force T 2 {k) to become trivial 
again. We have not been able to figure out any reasonable approximation which produces a 
non-trivial T 2 (k). 
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Anyway, in order to illustrate the kind of fixed point we are looking for, let us take 
another option which, unfortunately, is completely unrealistic. It consists of sending i0\ \K\ 
and CO3 (^77) to zero, keeping uj 2 (p?) as a mere constant (:= a 2 ). That presents the main 
advantage of producing decoupled equations for T\ and T 2 : 

T x = c (l + A), 

T 2 = ci(l + a 2 B). (C.9) 

where A := T\ X and B = T 2 X . Both are well defined provided c (l + A) and ci(l + a 2 B) 
are finite. We compute them multiplying above by 1/(E — jj -\- irj) and integrating. This 
produces: 

1 

c (l + A) 

C!(l + a 2 B) = ^— . (CIO) 

Tl -a> 2 l 

It is obvious that divergences are absorbed if Cq, Ci behave like A -1 and non-trivial results 
(Ti, T 2 7^ 0) require: 

1 MA 1 

+ 



_1_ 

CO 


-2b' 

1 


Cl 


-a 2 2 Q 



c 2vr 2 c[S(/i) ' 

- - _M^l + J_ (C11) 

Namely, c\ must be fine-tuned (to a non-trivial fixed point) as desired. Unfortunately, 
as mentioned before, the assumptions made for the a;, here are not realistic. 

Summarizing, we are rather pessimistic about the possibility that a non-trivial RG fixed 
point for both Cq and c± exists, which allows for partial wave mixing at leading order. 

D. Proof that no continuous solution of eq. (18) of |{tgj 
exists when R — ► 

Consider eq. (18) of @, 



-MVoRcot(^MVoR) = ^ + ^-^ tan(2^/^p + O ) , (D.l) 

with Vq < , R > 0, 0o G [— 7r/2,7r/2]. We are interested in whether continuous solutions 
^o = K)(-R) exist when i? — ► 0. Let us define: 



y := ^/-MV R > , x:=2y^^ + o . (D.2) 
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In terms of these variables we are interested in whether a continuous solution y = y(x) exists 
when x — > oo for the following equation, 

3 x — 0o /T ^ 

ycoty = --\ — tanx. (D-3) 

Deriving this equation once one obtains: 

sin 2y — 2y dy sin 2x + 2x — 20o 



2 sin 2 y dx 4 cos 2 x 



[DA) 



which proves that y(x) decreases when x increases for x large enough. The proof holds 
everywhere except for the points x = (n + l/2)7r, y = inn, n, m = 0, 1, 2... which we analyze 
in the following. 

When x approaches (n + l/2)ir for a given n, y must necessarily approach rrnr for some 
m in order for eq. ( p.4|) to have a solution. If we write: 

x = {n + —)tx + 5x , y = m7r + 5y , &r , <5y — > , (D.5) 

we have, for m^O, 

2?77 

5y = ^^5x + 0(5x 2 ) (D.6) 

n+ ^ - tr 

Z 7T 

Hence, eq. ( |D.4| ) admits a continuous solution near the point x = (n + l/2)7r provided that 
we choose m^O. Notice also that y keeps decreasing when x increases in the neighborhood 
of this point. 

Now, if we increase x from (n + l/2)7r to (n + 3/2)TT, y must decrease from ran to (m— 1)tt, 
if continuity is required. By iterating the argument, if we increase x till (n + m + l/2)7r, 
continuity requires y to decrease till 0. However for x = (n + m + l/2)7r + 5x (5x — > 0) and 
y = by — > 0, eq. ( p.4| ) does not have a solution anymore, since one obtains: 

1 + 0( V) = - ( " + m 2 + ,f -^ (l) (D.7) 

This implies, in particular, that the curves plotted in fig. 4 of [1^] cannot be continuously 
extended below R ~ 0.25 fm., R ~ 0.13 fm. and i? ~ 0.09 fm. respectively, as it shown in 
the figure below. 

In conclusion, no continuous solution y = y{x) of eq. ( p.4|) (and hence of eq. (p.l|) ) 
exists for x — > oo (R — > 0). If continuity is given up, an infinite number of solutions exist, 
none of them being compatible with a RG flow, at least in the standard sense. Note also 



that this situation is qualitatively different from a limit cycle behavior [28], which is realized 



for instance, in three body systems [29|. There the flows are oscillating and discontinuous 



but uniquely defined, no matter how large the cut-off is. 
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Figure 1: In the figure it is shown how the first three branches of the flow presented in Figure 4 



of | IE], behave as one approaches the relevant limit R — ► 0. 
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